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ABSTRACT 



We report on the detection of a transit of the super-Earth 55 Cnce with warm Spitzer in IRAC's 4.5/im band. Our 
MCMC analysis includes an extensive modeling of the systematic effects affecting warm Spitzer photometry, and yields 
a transit depth of 410 ± 63 ppm, which translates to a planetary radius of 2.08lo J^^ ^© measured in IRAC 4.5/im 
channel. A planetary mass of 7.8llQ g3 is derived from an extensive set of radial-velocity data, yielding a mean 
planetary density of 4.78li['2o gem"''. Thanks to the brightness of its host star {V = 6, K = A), 55 Cnce is a unique 
target for the thorough characterization of a super-Earth orbiting around a solar-type star. 

Key words, binaries: eclipsing - planetary systems - stars: individual: 55 Cnc - techniques: photometric 



1. Introduction 

Radial velocity (RV) , microlensing and transit surveys have 
revealed the existence in our Galaxy of a large population 
of planets with a mass of a few to ~20 Earth masses (Lovis 
et al. 2009; Sumi et al. 2010; Borucki et al. 2011). Based on 
their mass (or minimum mass for RV planets), these plan- 
ets are loosely classified as "super-Earths" {Mp < 10 M^) 
and "Neptunes" {Mp > lOM^). This classification is based 
on the theoretical limit for gravitational capture of H/He, 
'^lO (e.g., Rafikov 2006), and thus implicitly assumes 
that Neptunes are predominantly ice giants with a signifi- 
cant H/He envelope, and that most super-Earths are mas- 
sive terrestrial planets. Still, the diversity of this planetary 
population is probably much larger than sketched by this 
simple division, as we can expect from the stochastic nature 
of planetary formation. 

The first transit of one of these low-mass planets, 
G J 436 b, was detected in 2007 (Gillon et al. 2007). Thanks 
to its transiting nature, the actual mass {Mp — 23.2 ± 
0.8 Me) and radius {Rp = 4.22±0.10i?e) of G J 436 b could 
be accurately determined (Torres, 2007), indicating for this 
"hot Neptune" a mass, radius and density indeed very sim- 
ilar to the ice giant planets Uranus and Neptune. More 
recently, several other transiting low-mass planets were de- 
tected. While many more planet candidates detected by 
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* The photometric time series used in this work are 
available in electronic form at the CDS via anony- 
mous ftp to cdsar c.u-strasbg.fr (130.79.128.5) or via 



http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/, 



the Kepler mission are waiting for confirmation (Borucki 
et al. 2011), the first confirmed low-mass transiting planets 
already show a large diversity. Some of these planets, like 
HAT-P-llb (Bakos et al. 2010) and Kepler-4b (Borucki et 
al. 2010b), are similar to Neptune and GJ436b. Kepler- 
11c (Lissauer et al. 2011) seems to be a smaller version of 
Neptune, while HAT-P-26b (Hartman et al. 2010) has a 
much lower density (0.4 ± 0.10 gcm^'^ vs 1.64 gcm^^ for 
Neptune) that is consistent with a significantly larger H/He 
fraction. The super-Earths CoRoT-7b (Leger et al. 2009, 
Hatzes et al. 2010) and Kepler-lOb (Batalha et al. 2011) 
are probably massive rocky planets formed in the inner part 
of their protoplanetary disks. The super-Earth GJ 1214b 
(Charbonneau et al. 2009) is still mysterious in nature. 
Its large radius {Rp = 2.44 ± 0.21 i?e. Carter et al. 2011) 
suggests a significant gaseous envelope that could originate 
from the outgassing of the rocky /icy surface material of a 
terrestrial planet or that could be of primordial origin, mak- 
ing it a kind of "mini-Neptune" (Rogers & Seager 2010). 
Recent transit transmission spectrophotometric measure- 
ments for GJ 1214 b seem to rule out a cloud-free atmo- 
sphere composed primarily of hydrogen (Bean et al. 2010, 
Desert et al. 2011), but more atmospheric measurements 
are needed to determine the exact nature of its envelope. 
The case of GJ 1214 b shows nicely that understanding the 
true nature of a low-mass exoplanet could require not only 
precise measurements of its mass and radius, but also a 
study of its atmospheric properties. 

Among all the low-mass transiting exoplanets detected 
so far, only G J 1214 b and GJ436b, and to a lesser ex- 
tent HAT-P-llb and HAT-P-26b, orbit around stars small 
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enough and bright enough in the infrared to make pos- 
sible a thorough atmospheric characterization with exist- 
ing or future facihties hke JWST (e.g., Shabram et al. 
2011). Improving our understanding of the low-mass planet 
population orbiting around solar-type stars requires that 
such planets are detected in transit in front of much 
nearer/brighter host stars than the targets of surveys like 
CoRoT (Barge et al. 2008) or Kepler (Borucki et al. 2010a). 
This is the main goal of two ambitious space mission 
projects in development: PLATO (Catala et al. 2010) and 
TESS (Ricker et al. 2010). Still, another and more straight- 
forward possibility exists. Doppler surveys target bright 
nearby stars, and they have now detected enough nearby 
low-mass planets to make highly probable that a few of 
them transit their parent stars. This motivated us to search 
with Spitzer for transits of low-mass Doppler planets having 
the highest transit probability. In a previous paper (Gillon 
et al. 2010, hereafer GIO), we described the reasons that 
have led us to conclude that Spitzer and its Infra- Red Array 
Camera (IRAC, Fazio et al. 2004) were the best instrumen- 
tal choice for this transit search, and presented the results 
of our Spitzer cycle 5 program targeting HD 40307 b (Mayor 
et al. 2009). The rest of our program consist of a cycle 6 
DDT program (ID 60027) of 100 hours that targeted ten 
other low-mass planets. Spitzer''s cryogen was depleted at 
the end of cycle 5, and these observations were thus carried 
out in non-cryogenic mode ( "warm Spitzer" ) . 

The recent announcement of 55 Cnc e transits detection 
by the MOST satellite (Winn et al. 2011) motivated the 
publication of this paper. Our initial analysis of the warm 
Spitzer data taken in last January concluded a transit de- 
tection but also that several sources of instrumental effects 
needed to be fully characterized before securing the detec- 
tion. We only recently obtained a satisfactory instrumental 
model for warm Spitzer photometry, through a global anal- 
ysis of calibration data and of all the observations of our 
cycle 6 program (Gillon et al., in prep.). Once applied to 
our 55 Cnc data, this instrumental model leads not only to 
the firm detection of the transit of 55 Cnc e, but also to a 
precise determination of its transit parameters. 

Section 2 presents our derivation of the transit 
ephemeris from the published RVs. In Section 3, we present 
our data and their analysis that reveals the transiting na- 
ture of the planet. We discuss our transit detection and its 
implications in Section 4. 

2. Transit ephemeris estimation 

We performed a global analysis of all the available RVs 
for 55 Cnc to estimate the most reliable transit ephemeris 
for 55 Cnce. This analysis was done with the adaptative 
Markov Chain Monte-Carlo (MCMC) algorithm described 
in GIO. We assumed Keplerian orbits for the five planets 
orbiting 55 Cnc, after having checked that planet-planet 
interactions had negligible influence on our solutions, us- 
ing for this purpose the Systemic Console software package 
(Meschiari et al. 2009). Our analysis was based on the or- 
bital solution recently presented for 55 Cnc e by Dawson 
& Fabrycky (2010). As shown by these authors, the or- 
bital period value initially reported for this planet, 2.8 days 
(McArthur et al. 2004; Fischer et al. 2008), was an afias of 
the true period, 0.74 day. We verified this result by making 
two independent MCMC analyses of the RVs, one assuming 
P ~ 0.74 day and the other assuming P 2.8 days. 



Using the Bayesian Information Criterion (BIG; e.g. 
Carlin & Louis 2008) to estimate the marginal likelihood 
of both models, and assuming that these models have the 
same prior probability, we obtained an odds ratio of ~ 10^^ 
in favor of the P = 0.74 day model, indicating a decisive 
strength of evidence for this model (Jeffreys 1961). The 
best-fitting model obtained from this analysis was used to 
estimate the jitter noise in the RV datasets. 6.0 ms^^for 
Lick, 4.3 ms~^for Keck, 5.5 ms~^for HET and 15 ms~^for 
ELODIE were added in quadrature to the published error 
bars to derive the uncertainties on the physical parameters 
of 55 Cnce presented in Table 1. 

In addition to some basic parameters for the host star, 
the origin of the RVs used as input data, and a descrip- 
tion of our warm Spiteer observations. Table 1 provides the 
most relevant results of our MCMC analysis for 55 Cnce. 
The large transit probability, ^-^29%, and the very well con- 
strained transit ephemeris (Icr error < 1 hour in 2011) of 
this super-Earth {Mp = 7.8±O.6M0) made it an extremely 
interesting target for our transit search program. 



3. Data analysis 

3.1. Data description 

55 Cnc was observed by Spitzer on 6 January 2011 from 
9h41 to 14h39 UT. The data consist of 5240 sets of 64 in- 
dividual subarray images obtained by the IRAC detector 
at 4.5 /im with an integration time of 0.01s, and calibrated 
by the Spitzer pipeline version S18.18.0. They are avail- 
able on the Spitzer Heritage Archive databasqj under the 
form of 5240 Basic Calibrated Data (BCD) files. We first 
converted fluxes from the Spitzer units of specific intensity 
(MJy/sr) to photon counts, then performed aperture pho- 
tometry on each subarray image with the IRAF/DAOPHOlQ 
software (Stetson, 1987). We tested different aperture radii 
and background annuli, the best result being obtained with 
an aperture radius of 3 pixels and a background annulus 
extending from 11 to 15.5 pixels from the PSF center. The 
center of the PSF was measured by fltting a Gaussian pro- 
file to each image. We discarded the first ten minutes of 
data to allow the detector to stabilize. The x-y distribution 
of the measurements was then looked at, and we discarded 
the few measurements having a very different position than 
the bulk of the data. For each block of 64 subarray images, 
we then discarded measurements with discrepant values of 
flux, background, x and y positions using a a median clip- 
ping (5tT for the flux and lOcr for the other parameters), and 
the resulting values were averaged, the photometric error 
being taken as the error on the average flux measurement. 
At this stage, a 50cr slipping median clipping was used on 
the resulting light curve to discard totally discrepant fluxes. 

Figure 1 shows the resulting raw light curve, and the 
time-series for the background and the x and y positions. As 
can be seen in Fig. 1 and Fig. 2, the measured background 
showed an unusual evolution during the run. It remained 
stable during ^3.5 hrs, then it increased abruptly of a few 
%, and finally its scatter increased largely. Such a behavior 



^ http:/ /sha.ipac. caltech.edu/applications/Spitzer/SHA 
^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Star 


55 Cnc 


Distance d [parsec] 


12.34 ±0.12^ 


V magnitude 


5.96 ±0.01^ 


K magnitude 


4.02 ± 0.03=^ 


Spectral type 


KOV - G8V* 


Effective temperature Teff [K] 


5234 ± 30^ 


Surface gravity log g 


4.45 ± O.OS'^ 


Metallicity Fe/H [dex] 


+0.31 ± 0.04^ 


Mass M, [Mq] 


0.905 ± 0.015'' 


Radius R* [Rq] 


0.943 ± 0.010® 


IW (lata 




250 Lick* 




70 Keck" 




119 HET'^ 




48 ELODIE 


Planet (MCMC results) 


55 Cnc e 


Minimal Mass Mp sin i [M©] 


7.80 ± 0.56 


Expected Radius Rp [-Re]" 


1.3 - 5.7 


Expected Area ratio {Rp/Rt)^ [ppm] 


150 - 3000 


Equilibrium temperature Teg [K]* 


1958 ± 15 


Ttransit - 2450000 [HJD] 


5568.011 ±0.025 


ToocultaUon - 2450000 [HJD] 


5568.368 ± 0.030 


Orbital period P [d] 


0.7365437 ± 0.0000052 


Central transit duration Wh=o [min] 


98 ±2 


RV semi-amplitude K [ms"'^] 


5.93 ± 0.42 


Semi-major axis a [AU] 


0.01544 ± 0.00009 


Eccentricity e 




Argument of periastron oj [deg] 




Prior transit probability [%] 


28.9 ± 1.5 


Prior occultation probability [%] 


29.3 ± 1.8 


Warm Spitzer data 


Channel [/xm] 


4.5 


AOR= 


39524608 


Exposure time [s] 


0.01 


Nbcd''' 


5240 


Duration [hr] 


5 



Table 1. Basic data for the star 55 Cnc, relevant results 
of our MCMC analysis of the RVs, and description of the 
data (RVs + warm Spitzer observations) used in this work. 
^Van Lccuwcn (2007). ^Turon ct al. (1993). ^Sknitskic ct 
al. (2006). ''Fischer et al. (2008). ^Valenti & Fischer (2005). 
''von Braun et al. (2011). '^Mac Arthur et al. (2004). ^Naef 
et al. (2004). "Assuming Mp sini = Mp. The minimum and 
maximum values correspond, respectively, to a pure iron 
and a pure hydrogen planet (Scagcr ct al. 2007). ''Assuming 
a null albedo and a heat distribution factor /' = 1/4 
(Seager 2010). ''AOR = Astronomical Observation Request 
= Spitzer observing sequence. ''BCD = Basic Calibrated 
Data = block of 64 subarray exposures. 
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Fig. 1. Top left: raw light curve obtained for 55 Cnc. Top 
right: background time-series for this run. Bottom: time- 
series for the X (left) and y (right) positions of the stellar 
center. 
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5568 5568.05 
BJD - 2450000 



Fig. 2. Background time-series for the 55 Cnc data, after 
normalization and binning per 15 minutes intervals. 



is most probably of instrumental origin. We included this 
instrumental effect in our data modeling (see below). 



3.2. Modeling the systematics 

The IRAC 3.6 and 4.5 /xm detectors are composed of InSb 
arrays that show a strong intrapixel quantum efficiency 
(QE) variability, the QE being maximal in the middle of 
the pixel and decreasing towards the edges. The full-width 
at half maximum (FWHM) of the point-spread function 
(PSF) is ~1.7 pixels. This undersampling of the PSF com- 
bined with the QE intrapixel variability leads to a strong 
dcpcndancc of the measured stellar flux on the exact lo- 
cation of the PSF center in a pixel. As Spitzer's pointing 
wobbles with an amplitude of ~0.1 pixel and a period of 
-^Ih, this leads to a severe systematic effect in the pho- 
tometric time-series acquired at 3.6 and 4.5 ^m, known 
as the "pixel-phase" effect. This effect was already present 
in the cryogenic part of the Spitzer mission and is very 
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5568 5568.05 
BJD - a45000U 



Fig. 3. 55 Cnc light curve corrected for the "pixel-phase" 
effect using a 2""^ order position polynomium and binned 
to intervals of 5 min. 

well-documented (e.g., Knutson et al. 2008 and references 
therein). It is the main limit to the photometric precision 
of warm Spitzer (e.g., Ballard et al. 2010). From a com- 
prehensive analysis, the Spitzer engineering team identified 
recently the cause of the Spitzer pointing wobble as the 
thermal cycling of a heater used to keep a battery within 
its nominal temperature rang€0. After extensive testing and 
review, it was decided to reduce by a factor of two the ther- 
mal amplitude of the cycling while increasing its frequency 
to make it differ more from the typical frequency of plane- 
tary transits and occultations. Our data were obtained after 
this heater change. The correlation between the measured 
fluxes and the stellar image position is clearly noticeable 
in the raw light curve, the resulting periodic pattern hav- 
ing a typical period ~ 35 min, corresponding to a cycle of 
the heater after the engineering change. We modeled this 
"pixel-phase" effect with the following 2'"'-order x and y 
position polynomial: 

A{dx, dy) fli -|- a2dx + a^dx^ + a4dy + a^dy^ 

+aQdxdy, (1) 

where dx and dy are the distance of the PSF center to the 
center of the pixel. This model for the "pixel-phase" effect 
is quite classical in the exoplanet literature (e.g., Knutson 
et al. 2008, Desert et al. 2009). Correcting the light curve 
with the best-fit "pixel-phase" model lead to the light curve 
visible in Fig. 3. It shows a drop of brightness with an am- 
plitude compatible with a transit of 55 Cnce. It also shows 
some other low-amplitude flux modulations that are caused 
by other warm Spitzer systematic effects (see below). 

One could argue that the transit-like pattern could be 
caused by the imperfect correction of the "pixel-phase" ef- 
fect by the function shown in Eq. 1. This is very unlikely, 
as the duration of the transit-like structure does not corre- 
spond to the one of the wobbles of the stellar position on the 
chip. To discard firmly this possibility, we tried 3'''* and 4*''- 
order version of Eq. 1 that led to very similar light curves. 
We also corrected the "pixel-phase" effect by a totally dif- 
ferent method that relies only on the data themselves and 
not on any numerical function. We divided the pixel area 
sampled by the PSF center into 33 x 33 small boxes. If at 
least 5 subarray measurements felt into a given box, and if 
these measurements sampled at least 0.14 days (70% of the 
duration of the run) the corresponding measurements were 



divided by their mean value. If these two conditions were 
not met for a given box, its measurements were discarded. 
The reduction procedure was then identical to the one de- 
scribed above. The light curve obtained after this correc- 
tion by an "intrapixel flatfleld" was totally similar (pattern, 
scatter) to the one visible in Fig. 3. To assess the depen- 
dancy of the observed transit-like structure on the details 
of the reduction procedure several independent reductions 
of the data were performed by four of us (M. G., B.-O. D., 
D. D., P. C), all using different reduction and detrending 
procedures. We also tested performing photometry on the 
5240 images resulting from the averaging of the 64 subarray 
images of each BCD flle, using a median filter to reject out- 
lying pixels. Finally, we inspected the light curves obtained 
without background subtraction. In all cases, the obtained 
light curves were very similar to the one shown in Fig. 3, 
confirming the independence of the obtained photometry 
on the details of the reduction procedure. 

At this stage, we performed a thorough MCMC analysis 
of our photometry to deduce the transit detection signifi- 
cance, using as input data the raw light curve obtained 
with an aperture of 3 pixels. Our model assumed a mass 
of 0.905 ± O.OISMq for 55 Cnc (von Braun et al. 2011), 
and a circular orbit with P = 0.7365437 days for 55 Cnce 
(Sect. 2). We used the model of Mandel & Agol (2002) for 
the transit, in addition to the following model for the pho- 
tometric variations of instrumental and stellar origin: 



A{dx, dy, dt) 



ai + a2dx + a-^dx^ + a^dy + a^dy^ 
+aedxdy + a^dt 



1 sm 



dt — Og 



-l-aii logdt + ai2 logdi , 



(2) 



^ http:/ /ssc. spitzer. caltech.edu/warmmission/news/ 
21oct2010memo.pdf 



where dt is the time elapsed since 2455568.05 BJD, i.e. the 
time at which the background increases sharply (Fig. 1 & 
2). The an and ai2 terms were only applied for dt > 0. The 
six first terms of this equation correspond to the "pixel- 
phase" model (Eq. 1). The purpose of the linear term in dt 
is to model a possible smooth variation of the stellar bright- 
ness. The other terms result from our extensive analysis of 
our entire set of warm Spitzer data (Gillon et al., in prep.) 
and of available calibration data that lead us to conclude 
to a low-amplitude periodic variability of the effective gain 
of the detector, its typical period lying between 30 and 60 
minutes and its amplitude being in average of a few dozens 
of ppm. Considering the challenging photometric precision 
required by our program, it is very important to take it 
into account, justifying the sine term in Eq. 2. We are cur- 
rently working with the Spitzer engineering team to find 
the origin of this periodic variation. We also notice that a 
"background explosion" such as the one affecting the last 
part of our data is correlated to a sharp increase of the ef- 
fective gain of the detector that is very well modeled by the 
last two terms of Eq. 2. The MCMC uses the whole dataset 
to simultaneously fit for the transit model and the baseline 
function presented in Eq. 2. 
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3.3. MCMC analysis and model comparison 

The following parameters were jump parameter^ in our 
analysis: the planet/star area ratio (Rp/Rs)"^, the transit 
width (from first to last contact) W , the impact parame- 
ter b = acosi/i?*, and the time of minimum light Tq. We 
assumed a uniform prior distribution for these jump pa- 
rameters, but we imposed a Gaussian prior for the stellar 
radius i?* based on R^ = 0.943 ± O.OlOi?0 (Table 1). We 
assumed a quadratic limb-darkening law with coefficients 
ui = 0.0706 and U2 = 0.1471. These values were drawn 
from the theoretical tables of Claret & Bloemen (2011) for 
the IRAC 4.5 ^m bandpass and for Tag= 5250 K, logg =4.5 
and [Fe/H]=-|-0.3. The 12 coefficients of the baseline mod- 
els (Eq. 2) were determined by least-square minimization 
at each steps of the Markov chains (see GIO and references 
therein for details) . The correlated noise present in the LC 
was taken into account as described by GIO, i.e., a scaling 
factor f5r was determined from the standard deviation of 
the binned and unbinned residuals of a preliminary MCMC 
analysis, and it was applied to the error bars. Several bin- 
ning intervals ranging from 10 to 90 minutes were tried in 
preliminary short Markov Chains, and the maximal value 
for /3r, 1.35, was used in our analysis. 

We performed two new MCMC analyses, one with a 
transit of 55 Cnce, and one without. Figure 4 shows the 
resulting best-fit transit model and its residuals. The odds 
ratio (Eq. 2 alone) vs (Eq. 2 -I- transit) is ^ 10* in favor of 
the transit model. The transit of 55 Cnc e is thus firmly 
detected. The period of the sinusoid (oio) derived from 
the MCMC is 51 minutes, well decoupled from the tran- 
sit duration (96 minutes) and significantly longer than the 
pixel-phase timescale (35 minutes). Its amplitude is 115±27 
ppm. We show in Fig. 5 the different contributions of the 
spatially- and time-dependent terms of Eq. 2. This shows 
that the time-dependent terms are well decoupled from the 
transit pattern. Table 2 presents the resulting transit and 
physical parameters and la error limits derived for 55 Cnc e. 



iRp/R,f [ppm] 
b — a cos i/R, [R^,] 
Transit width W [d] 
To - 2450000 [BJD] 
Rp/ Rt 
a/R* 
Inclination i [deg] 
Radius Rp [R<s\ 
Mass Mp [Me] 
Mean density pp [gcm~ 



410 ± 63 



0.16: 
0.0665 



5568.0265: 
0.0202tg 



3.517; 
87.3 
2.08 
7.81 
4.78 



13 
10 

,0011 

,0019 

-0.0015 

-0.0010 

,0015 

,0016 

041 

040 

.7 

.1 

16 

17 

58 

53 

31 

20 



Table 2. Median and la limits of the posterior distribu- 
tions derived for 55 Cnc e from our MCMC analysis of our 
warm Spitzer photometry. The mass and mean density are 
derived from the parameters in Table 1. 



0.0005 


-0.0005 



556B 5568.05 
B,ID - 2450000 



Fig. 4. Top: 55 Cnc light curve divided by the best-fit base- 
line model, binned to intervals of 5 min, with the best- 
fit transit model overimposed. Bottom: residuals of the fit 
binned to intervals of 5 min. 



2006) to obtain an additional estimation of the residual 
correlated noise. We used for this purpose the lightcurve 
corrected from the systematic effects described in Eq. 2. 
The resulting parameters are in excellent agreement with 
the ones derived from the MCMC analysis (Table 2), while 
their error bars are significantly smaller. This result indi- 
cates that the error budget is dominated by the uncertain- 
ties on the parameters of the complex baseline model, and 
not by the residual correlated noise. 

To test the robustness of our transit detection and of the 
resulting transit parameters, we performed ^ 10 additional 
MCMC analyses as described above, each of them assuming 
a different set of time-dependent terms presented in Eq. 2. 
We used a binned lightcurve (per 30s) for the purpose of this 
comparison to speed up the analysis. Table 3 presents the 
baseline model, derived depth, BIC (Bayesian Information 
Criterion) and Bayes factor obtained for 4 of those MCMC 
analyses. 

While none of these models revealed to be better than 
our nominal model for representing our warm Spitzer data 
(Bayes factor between 10^ and 10^^), each of these mod- 
els lead to a decisive detection of the transit of 55 Cnc e 
(Bayes factor between lO^*' and 10*^°). For all these alterna- 
tives models, the deduced values for the transit parameters 
agreed well with the ones deduced in our nominal analysis, 
except for the transit depth when the jump is not included. 
Nevertheless, our Bayesian model comparison makes these 
alternative models > 4 x 10^ times less probable than our 
nominal model. Table 3 also illustrates how the jump and 
the sinusoidal variation terms improve the baseline model. 
We thus not only conclude to our firm detection of a tran- 
sit of 55 Cnc e, but also to the robustness of the deduced 
results shown in Table 2. 



We also conducted a residual permutation bootstrap 
analysis, known as the prayer bead method (Gillon et al. 

Jump parameters are the parameters that are randomly per- 
turbed at each step of the MCMC. 



4. Discussion and conclusions 

4.1. Photometric precision of warm Spitzer 

Because we have to take into account three different instru- 
mental effects in addition to a possible smooth variation of 
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0.9995 



0.05 
a 
-0.05 



5567.95 



5567.95 



,0.001 




5568 



5568 



5568.05 



5568.05 



5568.1 



5568.1 



Fig. 5. Spatially- and time-dependent terms of the model 
function used in Eq. 2. Model terms are described as follows 
: tr is the transit model, p is the pixel phase correction (a2, 
03, 04, 05 and oe in Eq. 2), s is the sinusoidal (ag, ag and 
oio) and j is the jump model (an and 012). The linear 
trend aj is not represented. 



the stellar flux, the complexity of our photometric baseline 
model is large (12 free parameters, Eq. 2). This illustrates 
well the challenge of ultra-precise time-series IR photom- 
etry, especially with a detector that is no longer actively 
cooled. By modeling this baseline in addition to the tran- 
sit in our MCMC analysis, we naturally take into account 
its uncertainties and their impact on the deduced transit 
parameters. Despite the complexity of the baseline model, 
we reach a very good precision on these transit parameters. 
This is due not only to the extensive characterization of 
the warm Spitzer instrumental effects performed by the ex- 
oplanet community and the Spitzer Science Center, but also 
by the extremely high-cadence made possible by the IRAC 
subarray mode. Indeed, we have here more than 5,000 pho- 
tometric measurements to constrain thoroughly the 12+4 
parameters of our global model. We show here that warm 
Spitzer can not only detect an eclipse of a few hundreds 
of ppm (Ballard et al. 2010), it can also measure its depth 
with a precision of ~ 60 ppm, leading to the conclusion 
that this space telescope has still an important role to play 
for the detection and characterization of transiting planets. 

4.2. Planetary radius 

Our MCMC results (see Table 2) yield a planetary radius 
of 2.08lo;Jg ^® measured in IRAC 4.5/Ltm channel. The 
error bars are determined from the posterior distribution 
function produced by the MCMC and includes the error on 
the stellar radius. The current planetary radius uncertainty 
is dominated by the error on the transit depth. 

On its side, the radius of the star itself is now extremely 
well constrained, thanks to recent interferometric observa- 
tions of 55 Cancri performed by van Braun et al. (2011) 
using the CHARA array. The resulting updated stellar ra- 
dius value (Table 1) now yields a negligible contribution 
from the stellar radius to the planetary size uncertainty. 

In the preprint version of this paper, we reported a plan- 
etary radius 30% larger than the one initially obtained by 



Winn et al. (2011) in the visible with MOST. After we sub- 
mitted our paper, a new version of the Winn et al. (2011) 
analysis was made available that yields good agreement 
with our results (at the la level). 

4.3. Composition of 55 Cnce 

We used the internal structure model described in Valencia 
et al. (2010) suitable for rocky and gaseous planets. We 
considered four different rocky compositions that span the 
possible range in radius. The upper bound for the radius is 
set by the lightest rocky composition, which is one where 
there is no iron. A planet with a radius larger than this up- 
per limit necessarily has volatiles. The lower bound for the 
radius is set by a pure iron composition. Both extreme com- 
positions are unlikely to exist given that 1) iron, magnesium 
and silicate have similar condensation temperatures, with 
the latter two making up most of the mantle of the Earth 
(i.e. Mgo.9Feo.1O-l- 8102) and 2) a pure iron composition is 
unlikely even with maximal collisional stripping (Marcus et 
al. 2010). The other two rocky compositions are an Earth- 
like one (33% iron core, 67% silicate mantle with 0.1 of iron 
by mol) and a "super-Mercury" (63% iron core, 37% silicate 
mantle no iron). We also consider volatile compositions in 
which we added different amounts of II2O or hydrogen and 
helium (H-He) at an equilibrium temperature of ^2000 K 
above an earth-like nucleus. 

The data obtained in this study for the radius 
(2.08tg;i^ ^e), and mass (7.81^55;^^ M®) place the com- 
position of 55 Cnc e intersecting the threshold line between 
planets that necessarily require volatiles (above the "no- 
iron" line), and the ones that may be rocky (below the 
"no-iron" line). However, most of the combinations of mass 
and radius lie above the upper limit of a rocky planet, re- 
quiring that 55 Cnc e have volatiles in its composition. We 
find that an envelope of a few parts in 10,000 of H-He or of 
order of ^10% water above an Earth-like core can fit the 
data well. In Fig. 6 we show the mass-radius relationships 
for the different compositions considered and the different 
known transiting super-Earths. Based on the same argu- 
ments proposed by Valencia et al. (2010) for CoRoT-7b, 
the timescale for evaporation of a H-He envelope would be 
too short (~ a few million years) for it to be considered 
as a plausible composition, whereas the timescale for water 
evaporation is of the same order of magnitude than that 
of the age of the system a few billion years). Thus, ac- 
cording to the Spitzer data analysed in this study, we favor 
a composition of an envelope of supercritical water above 
a solid, perhaps earth-like, nucleus. The exact amount of 
volatiles will depend on the composition of the solid nu- 
cleus, with a reasonable estimate around '^15%. However, 
a pure rocky composition cannot be ruled out, in which 
case the planet would be depleted in iron with respect to 
Earth. 

Similarly, the data for 55 Cnc e reported by Winn et al. 
2011 also lies at the threshold of these two types of planets, 
albeit with a denser composition. 

Figure 7 shows the density as a function of mass of sev- 
eral transiting super-Earths. While CoRoT-7b and Kepler- 
10 b have practically the same composition, 55 Cnc e, with 
its similar effective temperature and mass, has a much 
lighter composition. It lies between the high-density "super- 
Mercuries" and the volatile-rich planets Kepler-llb and 
G J 1214 b. Within a small range of masses, 4-9 M^, the 
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Fig. 6. Mass- Radius relationship for 55 Cnc e. We show four 
different rocky compositions: no iron, Earth-hke (33% iron 
core, 67% sihcate mantle with 0.1 of iron by mol), super- 
Mercury (63% iron core, 37% silicate mantle no iron), and 
a pure iron planet. We consider two types of volatiles com- 
positions: 0.1-0.01% of H-He (pink), and 10-50% of water 
(blue) above an Earth-like nucleus. We show our data for 
55 Cnce (red cross with label (1)), that reported by Winn 
et al. (2011) (blue and label (2)), and the data for the 
known transiting hot super-Earths, Kepler-lOb (K-lOb - 
data from Batalha et al. 2011), and CoRoT-7b (C-7b data 
from Bruntt et al. 2010, Hatzes et al. 2011). Uranus and 
Neptune are shown for reference. 

known transiting super-Earths span a relatively large vari- 
ety in compositions. 
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P 


p + t 


p + s + t 


P + j + t 


p + s + j + t 


[ppm] 


590 ± 72 


665 ± 70 


683 ± 87 


428 ± 62 


410 ± 63 


BIG 


804 


758 


758 


746 


732 


Bayes factor 




9.8 X 10^ 


9.8 X 10'' 


3.9 X 10^^ 


4.3 X lO^'^ 



Table 3. Transit depth, Bayesian Information Criterion (BIC) and Bayes factor from the MCMC obtained for 5 different 
model baselines. Model terms are described as follows : p is the pixel phase correction (02, as, 04, 05 and ae in Eq. 2), t 

is the time-dependent linear trend (07), s is the sinusoidal (ag, ag and 0,10) and j is the jump model (an and 012). The 
Bayes factor given in the table is relative to the p model. Our adopted model described in Eq. 2 is the rightmost one. 



